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1 Introduction. Some history 



The problem of the high-order behaviour of the perturbation-expansion coefficients in field 
theory calculations has aroused new interest during the past two years. Particular attention 
has been paid to power corrections to QCD predictions for hard scattering processes. One 
can point out two reasons for this growing interest: 

• The theoretical problem of how physical observables can be reconstructed from their 
(often divergent) power expansions. 

• The pragmatic need to assess the usefulness of performing the extensive evaluations 
of multi-loop Feynman diagrams in QCD. Much effort has been devoted to the com- 
putation of higher-order QCD perturbative corrections; in some cases third-order 
approximations are known, and we now seem to be at the border of what can be 
carried out analytically or numerically in high-order perturbative calculations. So, 
why next order? If the series is an asymptotic one, the next order may represent no 
improvement with respect to the lower-order result. On the contrary, at a certain 
order it will lead to deterioration. 

There are good reasons to believe that the vanishing convergence radius of perturba- 
tion expansions of physically relevant quantities is a general feature of quantum theory; 
it has been proved in most of the theories and models considered. Bender and Wu [|I| 
showed in 1971 that perturbation theory for one-dimensional quantum mechanics with a 
polynomial potential is divergent, and they also discussed the very-large-order behaviour 
of the perturbation coefficients. In 1976, Lipatov p| obtained the same results for massless 
renormalizable scalar field theories. 

Brezin, Le Guillou and Zinn- Justin 0, applying independently the same method to 
anharmonic oscillations in quantum mechanics, were able to rederive and to generalize the 
results obtained by Bender and Wu. 

In QED, after the pioneering work of Dyson a number of papers appeared discussing 
the analyticity properties of the Green functions (some of them are cited in refs. and 
0). The growth like n! (where n is the order of approximation) has two sources: 

1. The number of diagrams grows like n!, each diagram giving a contribution of the 
order of 1. 

2. There are types of diagrams for which the amplitude itself grows like n\ [^. 
Gross and Periwal @] proved in 1988 that perturbation theory for the bosonic string 

diverges for any value of the coupling constant and is not Borel summable. 

The situation in QCD is particularly complex, not only because the expansion coeffi- 
cients behave like n! and are of non-alternating sign, but also due to the strong dependence 
of a truncated series on the renormalization prescription. Of particular interest are situa- 
tions in which the kinematic regime allows one quantity, say M (momentum, rest mass), 
to have a large value. Then the matrix element of a QCD operator O containing quark 
fields is represented by means of the operator product expansion in inverse powers of M: 

O = C^{M/^i)OM + i^C2(M/^)02(/x) + 0(1/M2), (1) 
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Table 1: High-order behaviour of perturbation expansion coefficients 



where /i is the renormalization scale, Oi are local operators of the theory (ordered by their 
dimension) and Ci{M/fi) are the expansion coefficients, to be calculated in perturbation 
theory. (As O is independent of the renormalization prescription, the renormalization- 
prescription dependence of the quantities on the right-hand side must mutually cancel.) 

Relation (^ allows for the separation of short- and long-distance effects in the process. 
Perturbation theory does not allow (as will be discussed below) an unambiguous computa- 
tion of the coefficients Cj(M/ fi). This generates an ambiguity of the order of 1/M (see |Q|) 
in the determination of Ci{M/fi), which implies that one should not include the 0{1/M) 
term in (|l|) until Ci{M/fi) has been correctly computed. Leaving aside for a while the 
question of what a correct computation of Ci(M//i) means, we conclude that the problem 
of a perturbative determination of Ci(M//i) has to be solved before passing to, or simul- 
taneously with, higher terms in the expansion (|I]). In my talk I shall focus on the former 
subject, referring for the latter topic to the recent review talk of ref. and references 
therein. 

Table 1 gives a survey of the large-order behaviour of the expansion coefficients in some 
typical theories and models. It is intended for ffist information and should not be used for 
systematic analyses because some important conditions or restrictions are not mentioned. 
(Let me also point out that references in the table as well as in the talk as a whole are 
often made not to the original papers but rather to reviews or more recent papers. As a 
result, a number of relevant valuable papers are not quoted; I apologize to their authors.) 
For details the reader is referred to the original papers and to the anthology by Le Guillou 
and Zinn- Justin [|| (which contains a list of references up to 1990). 



2 



The series are divergent in all the cases indicated, with vanishing radius of convergence. 
In this connection, two remarks are in order. The first is that the estimates are based only 
on certain subclasses of higher-order diagrams which, in the case of QED (QCD), are 
obtained by inserting an arbitrary number of fermion loops into the photon (gluon) lines 
of the lowest-order radiative correction. It is not known which additional contributions 
come from other diagrams; whether they are negligible or give the same or even a greater 
contribution, or finally whether they cause cancellations in the original estimate. 

The second remark I want to make is as follows. In most of the items of Table 1, 
the coefficients exhibit the same, factorial, high-order growth (items 1 - 5, 9 and 10). It 
would however be misleading to conclude that all these theories face equal divergence and 
ambiguity problems, which could be treated in the same way. Knowledge of the high- 
order behaviour shows only one part of the problem of a perturbation expansion, the other 
ones being those of summability and of uniqueness of the summability prescription. Some 
details are discussed in sections 2 and 3. 

2 Useful facts on power series 

Certain important facts on power series are overlooked in physical considerations. It may 
therefore be useful to recall some of them here because spontaneous intuition is often 
misleading. 

1. The divergence of a perturbation expansion does not signal an inconsistency or 
ambiguity in the theory. (See an analogue in item 5.) 

The problem is not that of convergence or divergence, but whether the expansion 
uniquely determines the function or not. The method of Feynman diagrams allows one 
to find, at least in principle, all coefficients of the perturbation series, which may deter- 
mine the function uniquely even if the series is divergent, and may not do so even if it is 
convergent. This depends on additional conditions. 

2. The requirement of asymptoticity of a perturbation series, 

oo 

f{z) ~ ar.z\ (2) 

n=0 

is not a formal assumption. It has physical content. 

When a perturbation series is divergent, it is usually reinterpreted as an asymptotic 
series. This is a weaker assumption, but not a technical one. It is by no means trivial; (0) 
for 2; — i> physically means that there is a very smooth transition between the system with 
interaction and the system without it. For certain classes of observables the perturbation 
series is believed to be an asymptotic expansion. 

3. If f{z) is singular at the origin, its asymptotic expansion (|^) may still be a convergent 
series. 

Consider for example the function f{z) = g(z) + Ae~"^^ with A real and a positive, 
where g{z) is analytic at the origin: the asymptotic expansion of the singular function f{z) 
in the right complex half-plane is a convergent series. This convergent series is asymptotic 
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to many functions (most of which are singular at the origin), but its values converge only 
to one of them, g{z). They do so inside the Taylor circle, which extends to the nearest 
singularity of g{z). 

4. A very violent behaviour of the expansion coefficients a„ at n — > oo might make us 
expect that no function with the property f{z) ~ Z^J^o'^"-^" would exist. This fear is not 
justified; it was proved by Borel and Carleman (see for details) that there are analytic 
functions corresponding to arbitrary asymptotic power series. 

5. The Borel non-summability of a perturbation expansion alone does not signal an 
inconsistency or ambiguity in the theory. 

The Borel procedure is just one of many possible summation methods and need not be 
applicable always and everywhere. The problem is to find a method (not necessarily the 
Borel one) which is appropriate for the case considered. 

6. If the ttn behave very violently at n — ^ cx) (so that the Borel series J2'^^=o has 
zero convergence radius) one might expect that it would be sufficient to replace n! in the 
denominator by a sequence 6„ that grows faster than n\, in order to reach a more efficient 
suppression of the If this is done, the price to pay usually is that stronger conditions on 
analyticity will be required for the summation procedure to be unambiguous. Analyticity 
of the expanded function must be examined simultaneously with the asymptotic expansion, 
otherwise the same series can be summed to different functions. 

We can conclude (and will elaborate below) that uniqueness of a summation procedure 
(in other words, recoverability of a function by means of its asymptotic series, see |Q) re- 
quires a balance between the high-order behaviour of the series and the analyticity domain 
of f{z). A violent high-order behaviour can lead to a unique definition only if "enough 
analyticity" is available. 



3 Analyticity vs. high-order behaviour: 
balance for uniqueness 

How to deal with divergent series and under which conditions a power series can uniquely 
determine the expanded function are questions of fundamental importance in quantum the- 
ory. Power expansions are badly needed in physics, but additional conditions are required 
to ensure that they have precise meaning. These additional conditions should refiect some 
physical features of the system. 

The fact that the expansion coefficients, which grow asymptotically like n\, are of con- 
stant, non-alternating sign, is the origin of most problems connected with the nonunique- 
ness of perturbative expansions in QCD. 

Let us discuss a simple example to illustrate this crucial point. Consider a generic 
quantity D, calculated in perturbation theory with coupling z, 

oo 

D{z) = Y^a^z^. (3) 

n=0 
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This can be rewritten as 



D{z) = y anz'\lln\) dte-'r . (4) 
n=o 

If the series (j^) has a non- vanishing convergence radius r, the integration in (|^) can be 
exchanged with the sum. If, on the other hand, the convergence radius is zero, r = 0, 
we can give the series meaning by exchanging the order of integration and summation. In 
either case we obtain 



oo 



D{z)= I dte-'Y.(^n^^ = r dte-'B{zt), (5) 



^■ 



where B{zt) is the Borel transform of D{z). Taking a„ = n\ (finite-order coefficients are 
irrelevant for the character of singularities) we obtain 



oo 



D{z)= / dte"*- -. (6) 

JO 1 — 

This integral does not exist for z positive, nor is the Borel sum of such a series defined. 
The summation can be defined in many ways; there are infinitely many functions with the 
asymptotic expansion S^q^'-^" • 

Note that the non-uniqueness of the summation prescription is not a mathematical 
difficulty; it rather signals lack (or insufficient use) of physics in the theory. There are 
two kinds of conditions for a function f{z) to be uniquely determined by its asymptotic 
expansion X^tT^"^"- '^^^ function must 

i) have a sufficiently large analyticity domain K 

ii) satisfy upper bounds (uniform in z and A^) on the remainder f{z)—Y,o o-nz"' for each 
N above a certain value. When K has a small opening angle at the origin, the inequalities 
must be sufficiently restrictive in order to reach uniqueness. If the angle is large enough, 
the condition ii) may be weakened. 

Let us sketch how this works in the case of the Borel summation method. The series 
(0) is called Borel summable if 

a) its Borel transform, 

oo 

B{t) = J2 «nt7^!, (7) 

n=0) 

converges inside some circle, \ t \< 6 , 6 > 0; 

b) B{t) has an analytic continuation to a neighbourhood of the positive real semi-axis 
Re t > 0, and 

c) the integral 

g{z) = - e-'/^Bit)dt, (8) 
z Jo 

called the Borel sum, converges for some z ^ 0. 

Nevanlinna P2| gave the following criterion of Borel summability: 
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Table 2: Analyticity vs. high orders: a balance is needed for uniqueness 

Let f{z) be analytic in the domain Kirf) defined by the inequality Re^ > ^ (with rj 
positive), a disc of radius \r] bisected by the positive real semi-axis and tangent to the 
imaginary axis (see Fig. la), and let f{z) have the asymptotic expansion (^. If the 
remainder Rn{z) after subtracting A^ terms from f{z), 

N-l 

Rm{z) = f{z) - J2 (^nZ"", (9) 
n=0 

is bounded by the inequality 

\RNiz)\ < Aa"" N\\z\^ (10) 

uniformly for all z G K{r]) and all A^ above some value Aq; then the sum is determined 
uniquely and has the form (§) for all z G K{r]). 

For other types of regions, similar theorems hold with modified regions. Details are 



exposed, along with references, in the review p3[i ^ survey of typical cases is given in 
Table 2, together with conditions for a unique determination of f{z) from its asymptotic 
expansion. 



The function /i(n) used in Table 2 is defined as follows (see |^ for theorems relevant 
to the wedge-shaped analyticity region) 



fi{n) = / exp(-e*)rdt. (11) 
Jo 



Note the double exponential function in the integrand, which implies a (Inn)" behaviour 
of the function fi{n) at large n. 

The disc K{ri), the "drop" K{ri, p) and the wedge W are depicted in Figs. 1 a, b, and c 
respectively. The general rule for the use of Table 2 is as follows: To ascertain uniqueness 
in the determination of a function f{z) out of its asymptotic series (0), one has to check if 
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(b) 




Figure 1: Three summation methods (see Table 2) for three different analyticity and 
boundedness domains: (a) the disc K{r]), (b) the drop K{ri,p), and (c) the wedge W. 
Crucial is not the size of the domain, but the opening angle at the origin. 



the conditions in all columns in one row are satisfied for the case considered; otherwise the 
function is either overdetermined or not uniquely determined. In QCD, where the function 
is, according to ||2^, analytic only in the wedge W (see row 4 of Table 2), but the a„ 
behave like n\ at large orders (row 2 of column 2), the function is not uniquely determined 
by its perturbation expansion and the situation calls for additional (non-perturbative) 
information. 

A simultaneous use of Table 1 and Table 2 can tell us to what extent physical observables 
and Green's functions can be reconstructed from the asymptotic series (0) in a theory. 
Taking, for instance, the QCD large-order behaviour, A^N'^Nl, from Table 1 (item 9), 
we obtain in Table 2 the uniform bound in the 2nd row to be valid on the whole disc 
K{ri) of the complex coupling constant plane as the condition for uniqueness. Of course, 
we do not expect uniqueness in perturbative QCD and there is therefore no surprise that 
QCD violates the balance required in Table 2; indeed, the actual region of analyticity of 
two-point Green functions in QCD is much smaller than K{7]), having the form of a horn 
with zero opening angle at the origin see row 4 of Table 2 and Fig. Ic. 

Table 2 also shows that the large-order behaviour of the coefficients a„ is not the only 
criterion of the (Borel or some other) summability of the series (^. Even a series with 
a very tame behaviour of the coefficients a„ may not be Borel summable, in spite of the 
radical suppression of the a„ by the Borel factors n\. An example is discussed in P 



^This non-perturbative result (see 't Hooft ||2^) is obtained by combining analyticity and unitarity 
of two-point Green functions in the complex momentum squared plane with asymptotic freedom. Of 
interest is a remark by Moroz ||2^ that the use of the Callan-Symanzik equation is not crucial in 't Hooft 's 
argument, which works also if theoretical evidence for asymptotic freedom is replaced by experimental 
evidence (if available) for a logarithmic dependence of on a sufficiently high mass scale. 
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4 Generalized Borel transforms 



1. The functions Bp{t) and M{t) defined in Table 2 are generalizations of the Borel 
transform, which can be used in the various situations listed in Table 1 to reduce non- 
uniqueness, provided some additional information is available. More about the properties 
of Bp{t) and M{t) can be found in ^ ^ ^ and in references therein. 



2. I will now discuss another type of generalization of the notion of Borel transform, 
p8| , p9| , ^ , which makes use of specific structures of singularities that are typical for QED 
and QCD. 

Let me first make a general remark. Until now I have been mostly exposing math- 
ematical methods. Now I will pass to some practical aspects of my talk, assuming that 
two-point Green's functions have special singularities, the renormalons, in the Borel plane, 
a structure that is now almost universally adopted. These ideas are based on various math- 
ematical models developed in the late 70's and early 80's, but nowadays these features are 
often considered as true features of Nature. With this reservation I am passing to this 
subject. 

The electromagnetic current-current correlation function is a useful example to explain 
a typical structure of singularities in the complex coupling-constant plane. Denoting this 
function 11'^'^, 

W = ij d^xe-'''^(0| T{f{x)f{0)) |0) (12) 

= {g'^q' - q^q^m-q') (13) 

and taking R, the ratio of the total cross section for e+e^ hadrons to that for e"'"e~ — >■ 
muon pairs, which is related to its imaginary part: 

i?(s) = 127rlmn(s + i0+), (14) 

we introduce a modified quantity 11 defined as 

n(g^) = -An'Q' (^^^ n(g^) (15) 

(where = — g^), to avoid inessential logarithmic terms. The perturbation expansion of 
n in the powers of the coupling constant has the form 



^ n=0 



(16) 

The dependence of 11 on ^^(Q^) exhibits a complex structure of singularities in the coupling 



constant complex plane at and around the origin |£4[. Their nature can be conveniently 
displayed when studied in the Borel plane as singularities of the corresponding Borel trans- 
form. Perturbation theory suggests the following structure of the singularities of the Borel 
transform lEPI: 
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(1) Instanton-anti-instanton pairs |^ generate equidistant singularities along the 
positive real axis starting at t = 4, for t = 41,1 = 1,2,.... Balitsky |3l| calculated the 
behaviour of i?e+e-^hadrons Hear t = 4 and found the leading I — I singularity to be a 
branch point of strength ^(A^f — A^), where and iVf is the number of colours and of 
flavours, respectively. 

(2) Ultraviolet renormalons are generated by contributions behaving as q+i ~ {—bQ/l)^k\, 
for I = 1,2,..., leading to singularities located at t = — 2//6o on the negative real axis. 
Near the first of the points, I = 1, the singularity is (feo^ + 2)"^"'"''', where 7 is related to the 
anomalous dimension of local operators of dimension 6. 

(3) Infrared renormalons are generated by contributions behaving as q+i ~ {hQ/l)^k\, 
I = 2,3,..., leading to singularities located at t = I/Bq. Near the first of the points, the 
singularity behaves as {hot — 4)-i-2A/6o^ 

Brown, Yaffe and Zhai and Beneke |2^, 29, 30] use the information about the structure 
of the first infrared renormalon. Expanding Im 11 and 11 in powers of the coupling constant 
with the expansion coefficients a„ and c„, respectively, and defining their respective Borel 
transforms 



-r(n + i) 



(17) 



and 



ciz) = £0 + E 



ncr, 



-r(n + i) 

they obtain, by comparing the expansion coefficients [[ 
A and C: 

A{z) = sm{boz)C{z) , 

which turns out to be a consequence of renormalization-group invariance pO |. 
modified Borel transform J-'{z) by 



T(z) = f ^(^ + ^^) 



the following relation between 

(19) 
Defining a 

(20) 



(thereby accounting for the first infrared renormalon), the authors of ||29[ and |^ consider 
the case of a general beta function, which they choose in such a scheme that its inverse 
contains two terms: 

l/(3{g') = -l/{bo9') + X/{bog'). (21) 

They find that, for this form of I3{g'^), the relation (^) remains valid also for A{z) and 
C{z), the modified (according to (pO|)) Borel transforms of Imll and 11 respectively: 

A{z) = sm{boz) C{z) . (22) 

It was already pointed out that the concept of renormalon is at present applied to 
concrete physical situations. This poses a topical problem: To what extent are renormalons 
physical concepts and to what extent are they just artefacts of our, maybe inadequate, 
formalism. 

Table 3 gives a survey of singularities in the Borel plane. 
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I — I pairs UV renormalons IR renormalons 
Position t = 4,8,12,... -2//6o, / = 1,2, ... 2//6o, / = 2,3, ... 

Strength f{Nf-N) {bot + 1)-^+^ (foot - 2)-i-2A/feo 

of the first ^ A = 61/60, ^0 

singularity 

Table 3: Singularities in the Borel plane for SU(A^) QCD 



5 Remarkable phenomena in the Borel plane 

A look at Table 3 reveals that the positions of the singularities in the Borel plane as well 
as their strength depend on N and A^f; they will thus change as the latter are varied. This 
is relevant to our discussion because the position and strength of the singularities nearest 
to the origin of the Borel plane determine the large-order properties of the perturbation 
expansion. This phenomenon was recently discussed by Lovett-Turner and Maxwell p2[ . 
Here 1 shall repeat a part of the analysis made in ref. and summarize a discussion I 
had with the second author on this topic. 

1. The instanton-anti-instanton pairs are covered by the infrared renormalons if Af = 
(mod 3). This follows from the condition 4n = 2Z/6o, where 60 = (HA^ — 2Af)/6. 

2. The first instanton-anti-instanton singularity (t = 4) disappears if A^f > A^ and 
A'f = A^ (mod 6). 

3. Take A^ = 3 and A^f = 15 as a special case of item 1. Then 60 = | and the l-th 
infrared renormalon coincides with the l-th 1 — 1 pair. Since the first renormalon, / = 1, 
does not exist, the leading singularity is the first I — I pair. 

4. If A^f = 16, 60 = I and the first infrared renormalon is located at t = 24, coinciding 
with the 6th I — I pair. Items 3 and 4 are examples of situations in which, contrary to 
common opinion, instantons play an important role in large-order behaviour. Analogous 
situations occur for different values of A^; it generally holds that when the number of 
flavours is sufficiently high (approaching the fiavour-saturation value), instantons become 
more important than renormalons. 

5. If, on the other hand, the number of flavours decreases, the importance of renor- 
malons increases. For Af = 0, 60 = and the density of renormalons is the highest; for 
3 colours, there are 9 infrared renormalons below the first I — I pair, which covers with the 
10th of them. 

6. The case A^f = 15, A^ = 3 is very special also from the point of view of the strength 
of the infrared renormalon singularity, whose power is 44 in this case. Because of this, 
the nearest infrared singularity disappears, but also the nearest instanton-anti-instanton 
singularity disappears, because 15 = 3 (mod 6), as is required in item 2 of this list. 

7. It is worth mentioning that if another plane than the Borel one is chosen (see Table 
2, p 7^ 1), all the singularities in the Borel plane either shrink to the origin or run away 
towards infinity. 



10 



Besides these general properties of the singularities in the Borel plane, considerable 
progress has been made in recent years in the knowledge of the singularities in the Borel 
plane in special physical systems and special functions. In particular, much is known in the 
case of heavy-light and heavy-heavy quark systems. So, in the large-A^f approximation 
there is a finite set of renormalon singularities ||33| and a discrete infinity of the renor- 
malon poles. Also, in the case of heavy-light quark-antiquark systems the structure of 
singularities in the Borel plane is known. 

The effects of renormalons have been studied also in application to heavy-quark physics: 
see for renormalons in heavy-quark pole mass, |3^, |3^ for renormalons in inclusive 
heavy-quark decay rates and ||3^ for the case of exclusive decay rates. The connection 



between renormalons and power divergences in heavy-quark physics was studied in [40 



6 Resummation of renormalon chains 

The renewed interest in calculating higher-order perturbative corrections and in examining 
the high-order behaviour of perturbative series is intimately related to the investigation of 
renormalization scale and scheme dependence of a truncated series as well as to attempts 
to estimate its uncalculated remainder. In the past various criteria for finding a suitable 
renormalization prescription were proposed; they are based on an estimate of the size of 
the remainder. Examples are Stevenson's principle of minimal sensitivity ||3^, Grunberg's 
notion of effective charge |3^, and the BLM method of scale setting by Brodsky, 



Lepage and Mackenzie. The BLM method, which is based on an analogy with QED, was 
further developed recently and I shall briefly comment on it. 

The BLM prescription is a method of estimating higher-order perturbative corrections 
of a physical quantity, provided that the first approximation is known. It consists in the 
use of some "average virtuality" as scale in the running coupling. Instead of working with 
fixed scale, 

as{Q^) I d'kF{k,Q), (23) 



one averages over the logarithm of the gluon momentum: 

«s(gkM) / d'kFik, Q) ^ as{Q') I d'k (^1 - ^asiQ') In ^ j F(A;, Q) . (24) 

This replacement amounts to accounting for higher-order terms in powers of a.siQ'^) by 
making use of the renormalization-group evolution 

a.{-k') = am t f^^^y" i-H-emr-' (25) 



n=l 



and retaining only the first two terms in the sum. This approach was recently generalized 
| 37| , PI , by introducing the running coupling constant aglk"^) directly into the vertices of 
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Feynman diagrams, with k being the momentum "flowing" through the hne of the virtual 
gluon. This modification means replacement of (|^) by 



a,(g*2) J d^kF{k,Q) = J d^kas{-k^)F{k,Q) 



(26) 



with a. 



X 



4vr/(/?oln 



"a2 



Note, however, that the beta function P{as{Q^)) is ap- 



proximated by its first term only. 

This method has been applied in phenomenology to various physical observables, such 
as T decay hadronic width and heavy-quark pole mass ||4^, semileptonic B-meson decay 
and the Drell-Yan process ^5\. The method makes maximal use of the information 
contained in the one-loop perturbative corrections combined with the one-loop running 
of the effective coupling, thereby providing a natural extension of the BLM scale-fixing 
prescription. 

Ellis et al. use Fade approximants to develop another method of resumming the 
QCD perturbative series. The authors test their method on various known QCD results 
and find that it works very well. 



7 Concluding remarks. Criticism 

A typical feature of the present status of the QCD perturbative corrections is the trend to 
avoid explicit calculation of higher-order corrections and, instead, to improve the result by 
making use of some of the additional information we may dispose of. Such information may 
be, for instance, the renormalization group invariance (which allows us to introduce the 
running coupling constant instead of the fixed one), analyticity, or the structure of singu- 
larities in the Borel plane (where the regular location of singularities, to be approximated 
by poles, suggests the use of Fade approximants). 

As already mentioned, the notion of renormalons, originally introduced and used to 
investigate interesting mathematical models, is now widely considered to have concrete 
background in physical phenomena. This universal belief meets with criticism that argues 
that the singularities in the Borel plane are nothing but products of a special choice of 
renormalization prescription [^]. Methods generalizing the scale-setting procedure devel- 
oped by Brodsky, Lepage and Mackenzie meet with criticism based on the argument 



that the approach is not fully independent of the choice of renormalization prescription. 
Further research will clarify the issue. 

It seems that the present effort in further developing the idea of renormalon and the 
corresponding formalism will be helpful in finding a language appropriate for physical 
ideas in the non-perturbative sector. Generalizations of the scale-setting procedures be- 
come valuable by implementing new physical information without calculating higher-order 
perturbative corrections (which is not only cumbersome but also doubtful, due to the 
divergence of the series), by using some additional, perturbatively independent informa- 
tion. This idea is not new, appearing in theoretical physics whenever technical difficulties 
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force one to look for methods allowing the exploitation of all information on the system, 
including whatever is not adequately taken into account by the existing formalism. 
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